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The long time behavior of the reduced time evolution operator for unstable multilevel systems 
is studied based on the N-level Friedrichs model in the presence of a zero energy resonance. The 
latter means the divergence of the resolvent at zero energy. Resorting to the technique developed 
by Jensen and Kato [Duke Math. J. 46, 583 (1979)], the zero energy resonance of this model is 
characterized by the zero energy eigenstate that does not belong to the Hilbert space. It is then 
shown that for some kinds of the rational form factors the logarithmically slow decay proportional 
to (logt)~^ of the reduced time evolution operator can be realized. 

I. INTRODUCTION 

The exponential decay of unstable systems has been a well-known law since the early days of quantum theory. The 
quantum description of those systems, however, allows deviation from exponential decay both at shorter and longer 
timesi than those times over which the exponential decay law dominates^*^ The short time deviation was actually 
found in a quantum tunneling experiment^ while the long time deviation seems still not to have been detected in any 
quantum systemi^ The main cause that hinders the detection is considered as the smallness of the deviation at such 
long times ^ 

In a recent study, a method enhancing the long time deviation was proposedi^i The decay of the unstable systems is 
theoretically modeled in the time evolution of the survival probability of unstable initial state. The survival probability 
is just the probability of finding the initial state in the state at a later time t. Since it is rewritten in a Fourier integral 
of the spectral function, its behavior at long times is determined by that of the spectral function near the threshold 
of the energy continuumiSi^ The essential aspect of the method is then distorting the spectral function from the 
Breit-Wigner form and dislocating its peak toward the threshold energy. Mathematically, this causes a divergence 
of the spectral function, i.e., the resolvent at the threshold. Then, it is expected that the exponential decay period 
disappears and the survival probability at long times is increased. A similar idea was also considered in a related 
contexti^ In addition, in the analysis of the Friedrichs modelSiiSi that is often used for the study on the decays of 
the unstable systems, the survival probability at long times sometimes exhibits a power decay law slower than that 
in cases of no divergence 

These facts remind the author of the zero energy resonance proposed by Jensen and Katoii4 According to them, 
such zero energy singularities are classified by the zero energy eigenstates of the total Hamiltonian that either belong 
to or do not belong to the Hilbert space. The cases where such eigenstates exist are called the exceptional cases; 
otherwise they are referred to as the regular case. The result in Ref. ^4 is concerned with the three-dimensional 
system of the one particle in short-range potentials, and they proved that the time evolution operator asymptotically 
decreases as 0(t~^/^) for the exceptional cases, that is slower than 0(i~'^/^) for the regular case. However, to the 
author's knowledge, the zero energy resonance for the Friedrichs model seems not to have been examined in the 
previous studies including Refs. ^^^^^^^in spite of the wide applicability of the model to the various physical 
systems&iiJ^ii^ ~ 

In the present paper, we examine the zero energy singularities of the resolvent at the threshold energy for the 
Friedrichs model from the viewpoint of the zero energy resonance^ and clarify how the asymptotic behavior of the 
survival probability at long times is affected. The Friedrichs model2ii£ describes the system of the finite discrete levels 
coupled with the continuous spectrum, in which the former can be interpreted as the unstable excited levels of atoms 
and the latter as the environmental electromagnetic fields >ii*i^*i& We emphasize that the model is not restricted to 
the single level caso'''i^i^i^°i^bi^ii3,i5ii6 j^^^^ rather, the TV-level casefiSiiii2ii2*2aiSi In addition, we assume that the 
square modulus of the form factors vanishes at zero energy with an integer power ^^tiSi^ however it is treated without 
restriction to a specific form to some extent. Furthermore, since we only consider the initial state spanned by the 
discrete states, it is sufficient for us to see the reduced resolvent R{z) that is just the restriction of the resolvent to the 
subspace spanned by the discrete states. Then, the Fourier integral of R{z) that we call the reduced time evolution 
operator U (t) enables us to calculate the survival probability. In fact it is expressed by the square modulus of the 
expectation value of U (t) in a given initial state. We first study the zero energy eigenstates of the model which either 
belong to or do not belong to the Hilbert space. It is then possible to estimate correctly the asymptotic behavior of 
R{z) at small energies both in the regular case and the exceptional cases. The latter cases are examined in detail only 
for the first kind, where only the zero energy eigenstate not belonging to the Hilbert space exists. On the basis of 
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this analysis, we can derive the long-time asymptotic formula for U{t) in those cases. In particular, the logarithmic 
decay proportional to {\ogt)~^ of U{t) is shown to occur in the exceptional case of the first kind for our form factors, 
which is extremely slower than the power decays in the regular case and in the exceptional case for another type of 
form factor liiiiSii^ These results are shown in Theorems I VII . 1 1 and I VII . 2l 

The organization of the paper is as follows. We first explain in Sec. |n] the A^-level Fricdrichs model with an 
appropriate Hilbert space, and then in Sec. IIIII wc introduce the reduced resolvent R{z). Section Hvl is devoted to the 
identification of zero energy eigenstates in this model. It is then possible to obtain the asymptotic expansion of R{z) 
at small energies in Sec. where we examine the regular and the exceptional case of the first kind. By making sure 
of the relation between R{z) and U (t) in Sec. IVII the asymptotic formula for U{t) in the regular and the exceptional 
case of the first kind are derived in Sec. IV ill . Concluding remarks are given in Sec. IVIIII 



II. HILBERT SPACE AND THE A^-LEVEL FRIEDRICHS MODEL 



We shall use bracket notation; however it can be understood in a standard treatment based on functional analysis 
as in Refs. ^0^17 . The Hilbert space describing the unstable multilevel systems is here defined by 

H:^C^ ®L\{0,oo)). (1) 

A vector |c) £ is expressed by |c) = X]^=i Cnl"-), where |n)'s are the orthonormal basis of C^, so that {n\n') = (5„„', 
where Snn' is Kronccker's delta. L^((0, oo)) is the Hilbert space of the square-integrable complex function |/) of the 
variable uj defined on (0,cx)), i.e., 

/•OO 

\f) eL\{0,^))<^ / I/Mpdccxcx). (2) 
Jo 

In a standard notation using the (generalized) eigenstate \uj) of the multiplication operator by lo, \ f) is nothing more 
than 

/>OC 

I/) = / f{uj)\oj)diu, (3) 
Jo 

where = S{lli — uj') and 6{u! — oj') is Dirac's delta. Then, an arbitrary vector j^I') e H composed of |c) G 

and I/) G i^((0, oo)) is denoted by 

l*>:=k) + l/>, (4) 
and the inner product between any two vectors \'^) and |$) E H is defined by2^ 



($|vl-) := {d\c) + {g\f) = ^ Cc„ + / g*{uj)fioj)dco, (5) 

n=l -^0 

where (*) denotes the complex conjugate and \^) = \d) + \g) with \d) E and \g) £ i^((0,oo)). In particular, the 



associated norm of j^") is ||^|| := ^yJW\W), which is ensured to be finite for all \'^) E Ti. 

Let us now introduce the A^-level Friedrichs model for a description of the decay of the unstable multilevel systems. 
The Hamiltonian H of this model is defined by 

H := Ho + XV, (6) 

where Hq is the free part and V the interaction part of H, respectively, and A G M is the coupling constant. Hq is 
defined by 

N ^oo 

Hq -.^ y2 ^n\n) {n\ + u;\uj){uj\duj, (7) 

n=l 

where ujn E M with uji < lj2 < ■ ■ ■ < (^n, and its action is prescribed by i?o|^) = 'l2n=i '-^nCn\n) + i^\f) for any 
l^f) = |c) + I/) E D{Ho). D{Ho) is the domain of Hq defined by D{Ho) := {|*) G Tl! | J^^ \uj f {uj)\^ duj < oo }, and 
then the self-adjointness of Hq is guaranteed. The interaction part V is defined by 

v-=T. KMN(c.| + «„HMH]d^, (8) 
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where we assumed that \vn) G L^((0, oo))i24 We call the L^-functions Vn{uj) the form factors of the system under 

consideratfon. The action of V is then given by = X]ri=i(""l/)l'^) + ^n=i ^n\vn) for any \^) G Ti. Note that 

since D{V) = Ti and F is a bounded self-adjoint operator, H is self-adjoint with the domain D{H) = D{Ho)nD{V) = 
D{Ho). 

In the whole of the paper, we will restrict ourselves to the special kind of the form factor: Suppose that the product 
f*j(w)t;„(a;) between an arbitrary pair of v^{ijj) and is written in a rational function, i.e., it is expressed by 

V^{uj)Vr,{uj) = — , (9) 

where ■Kmn{'^) and Pmn{^) are the polynomials of the degree A'/„i„ and Nnm, respectively, and we assume that pmn{^) 
has no zeros in [0, oo). It is also assumed that A/,„„ -f 2 < and 7rmn(0) = 0. The former condition ensures that 
v^(uj)vn{(jj) is integrable in [0, oo) and lim^_>oo w^(w)i^n(^^) = 0, while the latter condition implies that the rational 
function u*„(ti;)w„(w) = 0(w) as w — > 4-0. The form factors with such properties are often found in actual systems 
involving the process of the spontaneous emission of photons from the hydrogen atom^^*^ and quantum dotsi^^ We do 
not treat the algebraic form factor that behaves as 0(a;^/^) as lj ^ +0 instead, associated with the photodetachment 
of electrons from the negative ionSiiSiiiS^ and the spontaneous emission from the atoms in the photonic crystalsjii 
however, the discussion developed in the following could be easily extended to such a case. 



III. REDUCED RESOLVENT FOR THE TV-LEVEL FRIEDRICHS MODEL 



In the following, we introduce the reduced resolvent that is simply the restriction of the resolvent of H to the N 
dimensional subspace C''^ © {0}. Since we only consider the initial state belonging to this subspace, this restriction is 
sufficient for our study. In a technical sense, this treatment corresponds to the appropriate choice of a weighted Sobolev 
spaceii^i2£ In the later sections, we do not distinguish the vector in from that in © {0}. After introducing the 
reduced resolvent, we see the existence of the boundary values of the reduced resolvent on the positive real line. The 
large-energy behavior of the reduced resolvent is also examined, which is necessary for a rigorous estimation of the 
long time behavior of the reduced time evolution operator. 



A. Reduced resolvent 



The resolvent of Hq and that of H are defined by Ro{z) = {Hq — z)~^ and R{z) = {H — z)^^, respectively, where 
we assume that z € <C\{(t{Ho) U (j{H)). (t{Hq) (or (j{H)) is the spectrum of Hq (or H), i.e., the set of the singular 
points of Ro{z) (or R(z)). Then, we have 

R{z)-Ro{z) = -Ro{z)VR{z) (10) 
= -Ro{z)VRo{z) + Ro{z)VRoiz)VR{z). (11) 

From Eq. (|10|l . one obtains the equation R{z) = (1 -I- Rq{z)V)~^ Ro{z) , which is the starting point of the asymptotic 
expansion of R{z) for the short-range potential systemsii^ On the other hand, we instead start from Eq. (|ll|l to 
obtain 

[Ho-z- VRoiz)V]R{z) = 1 - VRoiz). (12) 

This equation can be solved for our model if we confine ourselves to the state subspace © {0}iiS In fact, from the 
fact that (n\VRo{z)\n') = for any and |n') e © {0}, Eq. lEJ) reads 

N 

[(W„ - z)Snm. - X'^Snra{z)]Rm7i'{z) = 5nn' , (13) 

m— 1 

where S{z) and R{z) are the N x N matrix defined with the matrix components 

Smn{z) := {m\VRo{z)V\n) = H ^Mf)^^lMdu;, and R^n{z) (m|i?(z)|n). (14) 

Jo LO - Z 

We call S{z) and R{z) the self energy and the reduced resolvent, respectively. Note that S{z) can be analytically 
defined for all z € C\[0,oo). For a later convenience, we also introduce the matrix Kq and K(z) by 

Kornn '■= {m\HQ\n) = ujnSmn, and Kmniz) [A'o - A^S'(z)]to„, (15) 
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respectively. Then, Eq. (|13|l is equivalent to 

[K{z) - z]R{z) = 1, Vz e C\(a(iJo) U a{H)), (16) 

which implies that det[i4:(z)-z]det[^(z)] = 1, so that det[i4:(z)-z] ^ and det[^(z)] ^ for all z G C\{a{HQ)Ua{H)). 
Thus, the inverse of K{z) — z exists, and we have 

R{z) = [K{z) ^ z]-\ yz eC\{a{Ho)Ua{H)). (17) 



B. The boundary values of R{z) and its large energy behavior 

From the assumption on the form factors, every v^{u!)vn{uj) is continued to the whole complex plane as a mero- 
morphic function which wc merely denote as u^(z)u„(z). It may have a finite number of poles. Then, it follows from 
Lemma FA. II that S{z) can be reduced to the form 

Siz) = 5(0) + Aiz) - (log(-z))r(z), (18) 

where we choose arg(— z) = arg(z) — tt and < arg(z) < 27r. The matrix r(z) is defined with the components 

r,nn{z) ■■= V*^{z)Vn{z), (19) 

and satisfies r{z) ^ as z in C. S'(O) is the limit of S{z) as z — > in C\[0,cx)), which turns out to be unique. 
Indeed, as we see from the Appendix in Ref. 2^, S,nn{0) ~ /o°° "m(^)^n('^)/'^ d,^- A{z) is then defined through 
Eq. (|18|l and becomes a Hermitian matrix for real lo, whose components arc the rational functions of z without any 
singularity on [0,cx)). By definition, A{z) satisfies A{z) ^ as z ^ 0. One sees that the boundary values of S{z) at 
the half line (0, oo) exist and satisfjiiS 

lim S{u; ± ie) = D{uj) ± 7rir(w), (20) 

where 

D{uj) -.^ 5(0) + A{uj) - (logtj)r(w). (21) 
The matrix D{lu) is just of the components 

D^^(.) P r <(->"(-O rf^^ (22) 

Jo LO' -UJ 

where P denotes the principal value of the integral. Note that both D{lj) and r{uj) are Hermitian matrices and 

r{uj) > 0. 

In all the discussion developed in the following, wc assume that 

dct[/v±(w) - w] 7^ 0, Vw > 0, (23) 

where we introduced 

K^{lj) := lim K{Lj±ie) = Ko- X'^D{uj)t >^^mr{Lj), Vw > 0. (24) 

It is worth noting that condition H23() is equivalent to the requirement of no positive eigenvalues of Hq, whose 
eigenstates are normalizable. Indeed, if det[K^{Lu) — lu] ~ for some w > 0, there is a non-zero vector \rj) = 
'^n=i Vnln) £ such that [K^{uj) — uj]\ri) = 0. Since both D{bS) and P{yS) are Hermitian matrices, the latter 
equation implies that 



{^\[Ko-u:-\'D{i^)\\^)^^ and {vi\r (u^M 



N 2 
Vn{uj)r]n 



= 0. (25) 



Note that the latter relation means that r{uj)\ri) = because r{oj) > 0. Thus, Eq. (|25|l implies that r{uj)\rj) = and 
[Ko-X^D{u;)]\Tj)=Lu\7j),i.e., 

N N 

Vn{uj)rjn = 0, and ^[w,n(^mn - X"^ D„in{u!)]rin = W?7m, (26) 



for all m = 1, . . . , A^. This is merely the condition for the existence of a positive eigenvalue to of H 
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Lemma III.l Under the assumption i2Sp . it holds that R,^{uj) := lim^^-i-Q R{uj ± ie) exists for all u} > and 
Proof : Under the assumption (|23|l . [K^{uj) — uj]~^ exists. Then 



\\[K^{lu)-u;]-^ -R{u;±ie)\\ 
< II [K^{uj) - uj]-^\\ \\±i€ + X^Siu ± ie) - \'^D{uj) T \^mr{uj)\\ \\R(u ± ie)|| . 



(27) 



Note that for any nonzero \y) e C (7^ 0) there is a nonzero |.t) e C such that \y) = [K{uj ± ie) — u ^ ie]\x). We 



then obtain 



\\R{uJ±^e)\y)\\ 

\\y\\ 



< 



< 



II [A'±(cj) - uj]\x)\\ - II [±ie + X^S{lu ± ie) - X^D{uj) T X^7rir{uj)]\x) 
||[A-±M-c.]|| 



inf 

-II ±ie + X^S{uj ± ie) - X^D{uj) T X^TTir{uj)\\ 



which imphes that 



hm \\R(uj±ie)\\ < 



inf 



||if±H-. 



< 00, 



(28) 



(29) 



where the norm of an iV x iV matrix A is defined by ||A|| = sup|^^_^Q i^-^gc" ll^l^)ll/l|2^ll- Eq. H28|l . we used the fact 
that there is some eo > such that for any positive e < eo and for any non zero \x) G 



> \\[K^iu;)-u;]\x 
" " ~ \x)^o,\x)eC'^ \\x\\ 



X 



> \\±ie + X^S{uj ± ie) - X-'Diuj) T AVir(w)|| 
II [±ie + X^S{uj ± ie) - X^D{lu) T X^Trir{Lu)]\x)\\ 



> 



(30) 



where the assumption (|23|l is taken into account. Thus, by using Eq. (EH), Eq. ^ leads us to 

hm II [K^{uj) - Lu]-^ - R{lu ± ie)\\ ^ 0, 

which completes the proof of the lemma. D 

Lemma III. 2 ; Under the assumption k2'i]) . R^{uj) is r-times differentiate in oj E (0, 00), and it behaves as 

d'-R^iuj] 



(31) 



duj^ 



= 0{uj ' ) as — > 00. 



(32) 



Proof : We first show the statement for r = 0. From the assumption on the form factors and Lemma lA.ll one sees 
that 



lim D{uj) = and lim r{uj) = 0. 



(33) 



Since from the assumption H23|l K^{uj) — a; is invertible for all > 0, it holds that there is some positive tj > wjv 
such that for any u) > ui, 



lly|| 



< 



< 



\{Ko - Lo)\x)\\ - X^\\[D{uj) ± nir{u:)]\x)\\ 
^ 0{u^-') 



u) — ujn — A^ll-D(ti') ± T:ir{uj)\\ 



(34) 
(35) 
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where the last inequality is obtained as follows: we can choose some positive lu > lon such that for any uj > to 

||[A'0 -Hl^'>ll • r 1 .2ur^, ^ , ^, xn .2\\[Dioj) ± X^TTir(Lo)]\x)\\ 



where Eq. was used. Thus Eq. (|35|l reads just as Eq. H32|) does for r = 0. In the case of r > 1. we first note 
that from our assumptions on the form factors and Lemma I A . II again . A{uj) and r{uj), which arc connected through 
D{ui) = 5*0 + A{uj) — \ogLjr{uj), also satisfy 

as w — !■ oo, where we used the estimation that fj!^"^ = 0{ijJ^^^''). Thus, for r — 1, we have 

^I^p^ ^ R±{oj)^ + X^D{iu) ± X'^TTiriuj)] R^ioj) = 0{uj-^), (38) 
as u; — > (X), where Eq. ^62\ for r ~ was used. For r > 1, we obtain 

^^^ = E E >HWLi){n^?^H£i-[- + A^^H±AW(..)]U±(c.), (39) 

where cS^\{si\l^i) is an appropriate positive integer. Note that the symbol ( ') means that the summation over 
{■5i}i=i is taken under the condition that > 1 for all i and X]i=i = r. If r = 1, Eq. (|39|) reproduces Eq. H38|) 
with a*-^H{'Si}i=i) = 1- the general case, if Eq. H39|l holds for r ^ k, then its derivative is made up of a linear 
combination of 

1 n ^^(^)t^ + ^'^(^) ± AVzr(a;)] i?±(c.), (40) 



where J2i=i Si = k + 1 for 1 < j < k, and 



|ri^^(^)£- + ^'^(^) ± A27rzr(c.)] I i?±(^), (41) 

where X]i=i Si = k + 1 for 1 < j < fc. On the other hand, they are actually included in the right-hand side (rhs) of 
Eq. ^ for r = fc + 1. Thus Eq. ^ is vahd for all integer r > 1. Let us now evaluate the asymptotic behavior of 
(f (uj) / du!^ for large uj. One can see that the summand for j = r in Eq. (|39() . where all Si = 1, contributes 0{uj~'^~^) 
to R^{lo)/ duj^ , while the other summands for j < r specified by {si}l^^ contribute 0(a;^'"^^^^'*" (logw)*")) at most, 
where sq is a number of Si satisfying Si > 2 and never vanishes for j < r. Therefore, the summand dominating for 
large uj is that for j ~ r. Since we recursively show a^^\{si\\^i) = r!, which never vanishes, the statement is proved. 
□ 



IV. CLASSIFICATION OF THE ZERO-ENERGY SINGULARITY OF k^{uj) 

In order to prescribe the zero energy resonance in the TV-level Friedrichs model, we should identify the zero energy 
eigenstates in this model which either belong to or do not belong to 7i. In the case of the short-range potential 
systems^ this task needs some elaborate examination with an appropriately extended Hilbcrt space. On the other 
hand, in our case, it is rather easily performed, as is seen in the following. 

Let us first see whether the eigenvector !■(/;) e of — A^S'(O), belonging to the zero eigenvalue, can be actually 
extended to the eigenvector of H belonging to the zero eigenvalue of H . If Y^) = \t/j) + \f) G D{H) C 7i is a zero 
eigenvector of H, it should satisfy i/j^") = 0, or equivalentljiSi 

N 

uj^ijn + X{vn\f)^Ofoi-n^l,...,N, and tj/(a') + A ^ Vni'«(w) = 0. (42) 

n=l 
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The latter equation of Eq. H42() is immediately solved as 



/(^)^_A^ii^^i^, (43) 

which should be square integrable because we intend to find l^*) in Ti. If this is the case, u)f{uj) G I/^((0,oo)), i.e., 
|\[') € D{H) is ensured, and the substitution of Eq. (|43|l into (wn|/) is safely done. Then, we find that the former 
equation of Eq. H42|l is nothing more than 

(i^o-A25(0))|V') = A-(O)IV) =0, (44) 

where K{Q) := K^{Q) = Kq — X^S{0). However, it is noted that such an /(w) associated with lip) is not necessarily 
square integrable. Hence, we shall decompose the zero eigenspace of K{0), denoted hy M = {\tp) G K{0)\^p) = 0}, 
into two kinds of subspaces: Mi = (Mq © Ms)-^ and M2 = {\tp) £ M\f{uj) £ L'^{{0,ooj)}. Here Mq = M-^, and 
denotes the orthogonal complement of the subspace D. In short we have = Mq Mi M2. Then, as is expected 
from the definition, we have 

Ml C {|^) £ M\f{uj) i L2((0, 00))}. (45) 

Note that in general the subset on the rhs of the above is not a subspace. We call the zero energy resonance 
(or merely zero resonance) of H if Mi is not empty. We also introduce the projection operators Qq, Qi, and Q2, 
associated with Mq, Mi, and M2, respectively. What we next do is to introduce the terminology following the study 
of Jensen and Katopi^ 



Definition IV. 1 We call the system a regular case if it holds that ^ (j(K{0)), i.e., 

det[is:(0)] ^ 0. (46) 

In this case, is said to be a regular point for H . 

Definition IV. 2 We call the system the exceptional case if, instead of Eq. ()46|l . it holds that £ cr(/v (0)), i.e., 

det[K{0)] = 0. (47) 

In particular, if is a resonance but not an eigenvalue {Qi 7^ 0, (32 = 0), is said to be an exceptional point for H of 
the first kind. If is not a resonance, but an eigenvalue {Qi = 0, (52 ^ 0), is said to be an exceptional point of the 
second kind. If is both a resonance and an eigenvalue {Qi ^ 0, Q2 ^ 0), is said to be an exceptional point of the 
third kind. 



We here remark that in general a non-trivial solution of Eq. (|44f) does not exist, however we can find a special case 
where such a solution surely exists. Suppose that eigenvalues uj„ of Hq are positive, and all form factors w„(w) 
satisfying the assumption are linearly independent. Then increasing A gradually form to 00, we can find some 
critical values of A for which A'(0) has the zero eigenvalue. Let us denote the n-th eigenvalue of K{0) by k,i(0) where 
Ki(0) < ^2(0) < • • • < Kjv(O). Then, k„(0) turns out to satisfy the inequality 

w„ - AVjv(O) < Kn{0) <OJn- X^<Ji{Q), (48) 

where both of (Ti(0) and (7n{0) are positive constants and ensured not to vanishiSi Thus, for a sufficiently small |A| 
Kn(0) for each n > iV— Af+ + 1 should be positive, while for a sufficiently large |A| they should be negative. Furthermore, 
one easily sees that all k„(0) are continuous functions of A^. Therefore, we conclude from the intermediate value 
theorem that there is at least one critical value of A to make k„(0) = for each n > N — A^_|_ + I. We can actually 
find such special values of A in Fig. 1 depicted in Ref. The example mentioned here could be treated in a more 
general way with resort to the analytic Fredholm theoremSSiSi which tells us at most a finite number of the critical 
values exists. 

It is also worth remarking that the existence of the zero energy eigenstates that either belong to or not to the 
Hilbert space necessarily prescribes the small energy behavior of the form factors in the following way. Remember 
that under the assumption on the form factors, r{uj) defined by Eq. H19|l has an asymptotic form like 

TV 

r(t.) = 5^^"r„ + 0(c.^+i), (49) 

n=l 
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as w ^ +0. Then, if € A4i exists, it should satisfy 

(V'|A|V'>^0. (50) 

In fact, if (V'lAlV') = 0, we see that fiu) in Eq. (gSJ has to satisfy |/(cj)P = A2(-0|r(tj)|?/>)/tj2 = 0{l) as lj +0, 
however which concludes that /(w) is square integrable. This contradicts the assumption that \^^) g Mi. In order 
to make the condition H50|l be satisfied, at least Fi should not vanish identically. We can find such form factors in 
the physical systems for the spontaneous emission process of photons from the Hydrogen atomiS*S^ and the quantum 
dot^ On the other hand, the discussion mentioned here immediately implies the fact that if € M2 exists, this 
time it should satisfy 

(^|A|^)=0, (51) 

which just ensures the requirement that /(w) £ L^((0, 00)). However, note that Eq. H51|) does not imply Fi = {) 
identically and only requires A = on the subspace M2. 

A. The small-energy behavior of R^iio) in the regular case 

In this case, the same as in Eq. we can show that -R^(O) — {K{Q))^^. Furthermore, we can choose some 

positive 0^0 > such that 

\\{K{Q)y^\\\\[uj + \^D{oj)±X^^iF{Lo)-X^S{Q)]\\ < 1, (52) 
for all positive lu < uiq. Then, R^{uj) is expanded as a Neumann series, 

R^{uj) ^ {K{Q)[1- {K{0)y^[uj + X^D{oj)±X^TriF{oj) - X^S{Q)]]}-^ = lim Sn{oj), (53) 

where 

N 

Sn{uj) = ^{(X(0))-i[t^ + X'^D{uj) ± X\iF{cj) - X^S{0)]y{K{0))-\ (54) 

3=0 

for all positive ui < uiq with {(A'(0))^^[w + X^D{uj) ± X^niF{Lo) — A^5'(0)]}" = 1. Under our assumptions on the form 
factors, A(uj) defined in Eq. (|21|l is asymptotically expanded as 

A(co) = ^c^"A„ + 0(t^^+i), (55) 

n=l 

as w ^ 0. By using Eqs. H49|l and (|55|) . it also follows that 

D{uj) = 5(0) -tjlogwri +0(^2 log tj), (56) 

as w — > +0. Then, Eq. (|53|l tells us the dominant asymptotic behavior of R^iuj) becomes 

R^{uj) = (K(0))-i + 0(cj logtj), (57) 
as w +0, where (/C(0))~^ never vanishes in the regular case. 

B. The small-energy behavior of R{_z) in the exceptional case of the first kind 

In the exceptional case of the first kind, from the definition, Qi 7^ while Q2 = 0, so that Qo + Qi = 1- Then, 
i?(z) is divided into the following four terms, 

R{z) = QoR{z)Qi + QoR{z)Qo + QiR{z)Qo + QiR{z)Qi. (58) 

We now introduce the four matrices, 

Eki{z) = Qk[K{z) - z]Qi = Qk[Ko - z - X^SifS) + A{z) ~ {\ogz)F{z) + iT:F{z)\\Qu (59) 
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where k,l = 0, 1, and log(— z) — logz = —in is used. From the relation that [K{z) — z]R{z) = 1, they satisfy 

EkoQoRQi + EkiQiRQi = QkSku (60) 

for k,l = 0, 1. To solve the above equations we need to check whether En and Eqq are invcrtiblc in the subspaccs Mi 
and Mo, respectively. By using Eq. H49(l . Eii{z) is rewritten as 

Eniz) = \^z{\ogz)QiriQi - zQi - \^Qi{Aiz) - (logz)[r(z) - zA] + ^7^^(z)}Ql, (61) 

where QiK{0)Qi = is used. Note that A{z) — 0{z) and r{z) — zFi — 0{z'^) for our form factors, so that all terms 
excepting the first one of the rhs of Eq. H61|l arc of the order of 0{z). Furthermore, the exceptional case of the first 
kind imposes the fact that QiAQi 7^ [see Eq. (|5U|l ]. and QiAQi is positive definite in Mi and thus invertible in 
Ml. Hence, -Eii(w) is invertible for sufficiently small \z\ > 0, and the inverse can be expanded by the Neumann series 
as, 

00 

En\^) = E(^ii(~^))'t^1 (QiAQi)-^ (62) 

^-^ A^zlogz 

j=0 ^ 

^ {QiEiQi)-' +0{z-\logz)-^)^Oiz-\logz)-^), (63) 



A^zlogz 
for small \z\, where we define 

^11 W ■= iQiriQi)-'{zQi + A2Qi{A(z) - (logz)[r(z) - zA] + ^^^(z)}Ql}, (64) 

A^zlogz ^ 

which behaves as 1/ logz as z ^ 0. For £^00 we have, 

Eoo{z) = QoK{0)Qo - Qo[z + X''[A{z) - (log z)r(z) + inr{z)]]Qo, (65) 

where the first term of the above is invertible in A/q, and the last term vanishes as z — > 0. Hence, Eqq{z) is invertible 
in AIq for sufficiently small |z| > and is expanded as 

00 

Eoo\z)^-Y,{[QoK{Q)Qo]-'[z + X^[A{z)-{\ogz)^{z) + ^7^^{z)]]Qoy[QoK{Q)Qo]-' ^0{1), (66) 
j=o 

for small |z|. Furthermore, we obtain 

Ekiiz) = -X'QkiAiz) - (log z)r(z) +7rir(z)]Q, = O(zlogz), (67) 
for fc ^ ? as z ~> 0, because QoK{0)Qi = QiK{0)Qa = 0. Solving Eq. ^U^, we obtainS^ 

QoRQo = (Eoo-EoiE^^'Ew)-^ ^0{1), (68) 

QoRQi = —Eq^EqiQiRQi = —Eq^Eqi{Eii — EiqE^^qEqi)'^ ~ 0(1), (69) 

QiRQa = —QiRQiEioEf^Q = —{Ell — EiqEq^Eqi) ^ E^Eq^ = 0{1), (70) 

QiRQi = (Eii-EioE^^'Eoi)-' ^0{z-\logz)-^). (71) 

It is worth noting that since the relation [K{z) — z]R{z) = 1 is analytically continued to the second Riemann sheet 
through the cut [0, cx)), the above-mentioned results are also valid for such a continued region and the estimations 
obtained here can be applied without any corrections. 

When we consider the small energy behavior of R~{uj), it is convenient to expand £'11, differently from Eq. H61|l . 

as 

Eii{z) = A2z(logz-27^^)Ql^lQl 

-zQi - A2gi{A(z) - (logz - 2^^)[^(z) - zFi] - iTrr{z)}Qi. (72) 

All the above-obtained results are only changed by replacing the term logz with logz — 27ri. Then, we can obtain 
from Eq. 

E+{uj)= lim Eii{uj + ie) ^ X^\ogLj)Qir{uj)Qi - LuQi - X^Qi[A{uj) - iTrr{uj)]Qi, (73) 
while from Eq. ifT^ 

E^^{uj)= lim Eii{oj - ie) ^ X^{\oguj)Qir{uj)Qi - LuQi - X^Qi[A{uj) + nrr{uj)]Qi. (74) 
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C. The small-energy behavior of R{z) in the exceptional case of the second kind 

In the exceptional case of the second kind, it fonows that Qi ~ 0, Q2 0; ^i^d Qo + Q2 = 1- Let us consider 
the asymptotic behavior of the reduced resolvent at small energies, which is written in the following form, R{z) = 
1=0 1 QkR{z)Qi. We now introduce the four matrices again, 

Ekiiz)^Qk[Kiz)^z]Qi, (75) 

where A:, ^ = 0, 2. From the relation that [K{z) — z]R{z) ~ 1, they satisfy that 

EkoQoRQi + Ek2Q2RQi = QkSki, (76) 

for k,l = 0, 2. This time, E22 and Eqq are invertible in A/2 and A/q, respectively. In fact, from Eqs. (|49() and H55|l we 
have 

E22{z) = ~zQ2[l + \^Ai)Q2 - \^Q2[A{z) ~ zAi - {\ogz)r[z) + iTTr{z)]Q2, (77) 

where Q2K{0)Q2 = was used. Note that since A(z) — zAi = 0{z^) and Q2AQ2 = [see Eq. H51|) ]. the second 
term of the rhs of Eq. (|77|l is of the order of 0(z^log2;). Furthermore, since Q2A1Q2 > from (52AQ2 = and 
Lemma [a. 21 Q2(l + ^^Ai)Q2 > and invertible in Af2. These facts bring us the fact that E22{z) is invertible in M2 
for sufficiently small z > 0, that is 



00 -. 

E22{^) = -Y.^-E22{z)y-[Q2{l + \^A,)Q2]-' (78) 
= i[Q2(l + A2Ai)Q2]"' + 0(logz) = 0(z-i), (79) 



where 



E22{z) := i[Q2(l + \^A^)Q2]-^\^Q2[A{z) - zA^ - (log z)r(z) + i^r[z)]Q2. (80) 



For £'00, we next have 

EM = QoA'(0)Qo -Qo[z + A2[^(z) - (log2)r(z) + mr{z)]]QQ, (81) 

where the first term of the above is invertible in Afo, and the last term vanishes as \z\ — > 0. Hence, Eqq{z) is invertible 
in Mq for sufficiently small \z\ > 0, and the inverse is obtained as a Neumann series. On the other hand. E20 and i?02 
behave as 

Eki{z) = Qk[-z[l + A'^i) - A2[A(z) - zA^ - [\ogz)r{z) + ^ir{z)]]Qi = 0[z), (82) 

for small \z\ where k ^ I. Solving Eqs. H76() as in Eqs. H68() to (|71|l . one sees that QqRQq — 0(1), for k,l — 0,2, 
except 

Q2RQ2 = {E22 - E2oE^o^Eo2)-' = Oiz-'), (83) 
as 2: 0. In particular, the last equation is expanded as, 

Q2R{z)Q2 = £[^2^ £^20^00' ^02] '^2^ = -^[Q2(l + >?Ai)Q2]-^ + O(logz), (84) 

for small |z|, where we used Eq. (|79|l . 

We now remark that the zero energy eigenspace of H denoted by Ao is completely characterized by A/2. That is, 
there is a bijection from Af2 ffi {0} to Aq. From the discussion concerning Eqs. (|42|1 . H43|l . and (|44|1 . for any \^) £ Aio, 
there is a vector £ AI2 ffi {0} such that 

I*) = |^> - A r ^ _ XR„iO)VM, (85) 

Jo ^ 

where V is restricted to ® {0} and Ro{0) is the (unbounded) multiphcation operator of 1/uj in i^((0, 00)). Then 
we see V\ip) G D{Ro{0)) because {ip) € Af2 © {0}. Thus 1 - XRo{0)V is well defined as an operator from Ar2 © {0} 
to H. Now, Eq. H85|) tells us that 1 — XRo{0)V is a surjection from M2 © {0} to Aq. On the other hand, for any 
I*) e Afo, if I*) = 0, i.e., = (*|*), Eq. ^ implies that = (*|*) > Therefore, 1 - \Ra{0)V is also an 

injection from Af2 © {0} to Afo, and the proof is completed. 
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D. The small-energy behavior of R^(uj) in the exceptional case of the third kind 

In the exceptional case of the third kind, from the definition, Qi ^ 0, Q2 ^ 0, and Qq + Qi + Q2 — 1- The reduced 
resolvent is written in the form, R^{uj) = QkR^{uj)Qi. This time, we need nine matrices, 



for k, 1 — 0,1, 2. From the relation that [K^{uj) — ijj]R^{uj) = 1, they satisfy that 



(86) 



(87) 



for k,l = 0,1,2. The asymptotic behaviors of £'^(0;) are essentially examined in the preceding subsections, except 
for Efi2{u}) and £'^j^(ijj). Then, i?j^(ci;) becomes 

E^^ioj) = -Lu\^QiAiQ2 

-\^Qi[A{uj) - ujAi - (logtj)[r(cj) -cjA] ±7ri[r(cj) -wA]]Q2 
= -uj\^QiAiQ2 + 0{uj'^\oguj), 



(89) 



where QiK{0)Q2 = 0, Q1Q2 = 0, and AQ2 = are used. The last relation follows from the fact that (52A(92 — 
and A > 0. In addition, since Qir{iv)Q2 — 0{lj'^), we see that Q1A1Q2 = /g Qir{u))Q2to~'^dijj. By the same way, 
we also see that 



Ef^{uj) = -u\^QiAiQ2 + 0{uj^\oguj). 
To solve Eqs. (|S7|) . let us now put the N x N matrix E as 



and partition it into A 



, B 



Eqq Eqi 

EiQ Ell 

formula again, E^^ (= i?) is expressed as^ 



E : 



E02 

El2 



Eqo Eqi E02 
Eio Ell E12 
E20 E21 E22 



(90) 



(91) 



C = [ £'20 E21 ] , D = [ i?22 ] • Then, from the inverse matrix 



E-i 



[A-BD-iq-i -A-iB[D - CA-^B]-! 
[D - CA-iB]-iCA-i [D-CA-iB]-! 



(92) 



The validities of A ^ and D ^ are already ensured in the exceptional cases of the first and second kinds, respectively. 
Then, one sees that since A^^ = 0{uj-^ (loguj)-^), B = 0{bj), C ^ 0{uj), and D^^ = 0{uj-^), it holds that 
A-iBD-iC = 0((loga;)-i). Thus, [A - BD-iq-i exists for small w and [A- BD^iq-i = 0(w~i(logw)-i). We 
also show that [D — CA~^B]~^ exists for small uj and [D — CA~^B]~^ = 0{uj^^). To obtain the asymptotic forms 
of the matrix components of E^^ explicitly, some redundant calculation is required; however, it could be achieved by 
a manner as similar to that used in the preceding subsections. 

V. ASYMPTOTIC EXPANSION OF THE REDUCED RESOLVENT AT SMALL z 

We examine the small-energy behavior of the reduced resolvent only for the regular case and the exceptional case of 
the first kind. This analysis is crucial for determining the asymptotic behavior of the reduced time evolution operator 
at long times. 



A. The regular case 

Here, we introduce A{uj) :— lu/X^ + A{uj) and suppose that A{uj) and r{uj) behave as 



ria+N 



AiLo) 



A2 



(93) 
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as a; ^ 0, respectively, that is, An = for all n < Ua and = for all n < n^, while An^ ^ and ^ 0. Then, 
we obtain 

^uj + D{uj) = S{0) + A{uj) - {\ogLj)F{uj) = 5(0) - w"" log wr„, + w""i„„ + 0(a;/iH), (94) 
A 



as ^ +0, where 



uj"-^ logw (rtb < Ua) 

LO"-" (Ub > Ua) 



(95) 



It is important to note that the values of two parameters Ua and Tif, are not determined independently. We shall 
here consider rib as a controllable one. Wc first note that if > 2 then ria = 1 should be concluded, because from 
Lemma Fa. 21 we have Ai > 0, so that Ai = 1 / + Ai > holds. Therefore, the conditions jif, < Ua and rib > na can 
be realized only in the situations 

rif, = 1 and ria > I, and rib > 2 and ria = 1, (96) 

respectively. 

Lemma V.l ; Assume that is a regular point for H . Then the r-th derivative of R^{uj) asymptotically behaves 



0(1) 



(r = 0) 



0(1) 



0{lu^-'- (log Lup^-"-^) (r > 1) 



(r = 0) 
(1 < r < Ub) , (97) 



0(w"''-''(logw)^("''-'')) {nb<r< 2nb) 



for = 1, or Ub > 2, respectively, as lu 0, where = max{a;,0} and 9(x) ~ 1 for x > Q or Q for x < Q. In 
addition, the r-th derivative of R^{lj) is approximated by that of a finite series 



{K{0))-' + (A'(0))-i [-L^"^(loga;)A2r„, +a."»A2i„„ ± AViw"^r„,J (A'(0))-\ 
that is, it is shown that 

doj'-X ^ ' 

-(A(0))-i - (A(0))-i [-c^"'(logc^)A2r„, +c^"»A2i„„ ± AVic^"'^,] (A'(0))-i 

= 0(tj2-'-(logtj)i+«(2"'-)) (r>0) 

r ©(c.I^-'-]^) (0<r<nb) 

\ 0(w"''+i-'^(logw)^("^+i-'')) (jib + 1 < r < 2nf,) ' 

for Hb = I, or Ub > 2, respectively, as to 0. Here, na is restricted to the condition ^96]). 

Proof : The left-hand side (Ihs) of Eq. H99|l is written as follows: 

-(A(0))-i - {K{0))-^ [-c^"''(logu;)A2r„, +u;"»A2i„„ ± A^WaJ (A(0))-i 



< 



(A(0))-i - (A(0))-i -u;"^(logc.)A2r„, + lo^-X^A^^ ± Wilo^^F^, (^(0))-^ 



(98) 



(99) 



(100) 
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where Sn{i^) is defined by Eq. (|54|l . When r = 0, the first term on the rhs of the above is estimated from the special 
case of the below for = 1, 



< 



\\lu + \^d{lo) ± \^Tiir{u) - a2s'(o)||^+i||(a'(o))-i||^+2 



1 - II (i^ (0))-i II II [w + A2D(w) ± A%ir(w) - A25(0)] II 
as w — > 0. When r > 1, instead we have 



o{h{ujr+'), (101) 



< 



(102) 



where Eq. 1)39(1 is used, and here > 1 and ~ r should be satisfied. Note that the first term on the rhs of 

Eq. p02|l appears only for r > 2, which is estimated in the following. In the following estimations, we temporarily 
forget the restriction H96|l and consider the two general cases: rib < na and Ub > Ua- In the case of Ub < Ua, we can 
obtain for r > 2 



j 



<E E '«^''^({^^K=i)||^^Hir o(^^"^-'-)rio((iog 

= 0(602""-'- (log a')''("'' + l-'')+^(2n, + l-r)^^ 

as w ^ 0. For Ua < rib, 



E E «*^^({^oLi)i?^Hn 



[lo + \^D{uj) ± \^TTir{uj)\ R^{uj) 



(2 < r < ria + Ub - 1) 

(9(^no+«,-rQQg^-)e(n„+«,-r)) (n„ + Tlfc < r < 2nb) 



as a; ^ 0. We here used that 



^[tj + \^D(uj) ± A27rir(tj) - \^sm 



= 0(^0""-'' (logw)'''""-'')), o: 
for rib < ria, or ria < rib, respectively, as ^ 0. Eq. (1105(1 follows from 



C)(t^K-'-]+) (0<r<7i6) 
©(^"''-'■(logtj)''^"'-'')) (r > rib) 



duj^ duj^ duj^ 

as 0. Incorporating Eqs. ((103(1 . ((104(1 . and Eq. II105I) . with 

d'^^±(a;) 



< 



E E >^(MLi)<|^^Mn^[^ + A'^(c.)±A'7rzr(c.)] [i?±(c.) 



R'^iuj)— \uj + \^D{lj) ± X^-Kirluj)] R^iuj) 



(103) 



(104) 



(105) 



(106) 



(107) 
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we have 



0(1) 



(r = 0) 



0(cj"''~'"(logw)''("''-'"')) (r > 1) 



0(1) (r = 0) 

0(wK"'-]+) (l<r<?ib) , 

0(w"''-'~(logw)''("''-'^)) K < r < 2??,f,) 



(108) 



for rih < ria, or Ua < nt, respectively, as a; — > 0. Then, the first part of the statement can be shown under the 
restriction (|96|1 . Let us next examine the second term on the rhs of Eq. (|l()2|l . which reads for r > 1, 



< 2 



duj 



+ \^D{lj) ± \^T:ir{Lj)\ 



0(w2"^-'^(logw)i+''("''"'^'), 



(9(^"„ + [n„-r]+-) {r<nb) 



(109) 



for rif, < Ua or ni, > respectively, as w — > 0. We here used Eq. H101() with = 0. Therefore, substituting Eqs. 
(ESIl, iP^ . and flTHH) into Eq. fllT^ . one has for r > 1, 



0(w2".-r(iQg^)l+eK + l-r)-)^ 



(j' < «a + - 1) 



(110) 



for rib < JT-a or rib > Ua, respectively, as w — *■ 0. Note that this estimation is also valid for r ~ because it reproduces 
Eq. ((TUT|l for N ^1. 

Let us now evaluate the last term in Eq. For r > 0, we have 

-(A'(0))-i - (X(0))-i [-L^"Hloga;)A2r„, +tu"»A2i„„ ± A27riu;"^r„,] {K{0))-^ 

<||TO)-i' 



duj' 



-0(^ (logo;) ^' \ 0(,,».+i-(logc.)''("^+i-'-)) (r>7ib + l) ' 

for rib < 'T-a or Ub > ria respectively, as w ^ for any r > 0. We here used that for r > 0, 

d-(r(o;)-o."^Aj ^ o(c.K+i-^l"), ^''aog^)(^H-^"-^^J ^ 0(,,".+i-(logc.)«(-+i-)) 
0, and so forth. Thus, setting Eqs. (|11U|1 and (|lll|l into Eq. (|10U|) . we conclude that 



(111) 



(112) 



as w 



(jr(o))-i - (/f(o))-i -c^"Hiogc^)A2r„, +c^""A2a„^ ± AVic^"'r„, (if(o))-i 



0(tj"f+i-'-(logw)^("''+i"'')) (r>0,?i,,>2) ' 

0(a;K+i-'^l^) (0<r <nb) 

(9(^"i+i-r(-iQg^-)e(n,+i-r)) (n^ + 1 < r < 2nt) 



(113) 
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for ri}, < Ua or ni, > Ua, respectively, as w — > 0. By taking into account the restriction H96() . we can show the last part 
of the lemma. D 

To estimate the long time behavior of the reduced time evolution operator, the above-mentioned lemma seems 
not precisely appropriate because the reduced time evolution operator is obtained from the Fourier transform of the 
imaginary part of i?+(a>), not from R'^{uj) itself, which is explained in the next section. Hence, the following lemma 
is more appropriate for our purpose. 

Lemma V.2 ; Assume thatO is a regular point for H . Then ther-th derivative oflmR'^iLu) :~ [R^ [lj) — {lo)) / 2i 
is approximated by that of 



in the sense that for < r < ni, + I the remainder is estimated as, 



(114) 



dw 



(115) 



for rif, = 1, or ujj > 2, respectively, as uj —^ 0. For r = + 2, the estimation is replaced by 0{uj ^) for ujj ~ \, 
Oiloguj) for Uf, = 2, or 0(1) for Uf, > 3, respectively, as cu 0. 



Proof: Since Im^+(tj) = X'^ttR+ {uj)r{uj)R- (cu), one has 

%{lmR+{co)^{K{0)r'X\to-^ir„,+cor„,+,+co'-r,,,+2)iK{0)r'} 



< A^TT 



J2 [Pi (s, t, u) + F2 (s, t, u) + Fsis, t, u)], 



(116) 
(117) 



s > 0,t > 0,u > 0, 
{s + t + u — r) 



with 



Fi{s,t,u) = Cstu 



d^R+iu 



F2{s,t,u)=C. 



duo 
d'R+{Lo) 



d*r{uj) 



d" 
dui 



-{R-{cj)-{Km-'} 



(118) 



duj'' 
d"(A^(0))-i 



dio 



j{r{uj)-uj"^{r^,+ujr^,+i + uj^rr,,+2)} 



duj-^ 



F'i{s,t,u) = Cstu 

d"(isr(0))-i 



£^{i?-(.)-(/.(0))-} 



du)" 



(119) 



(120) 



where CstuS are appropriate constants. For 1 < r < nt + 1, the summation of the first summand in Eq. H117|) can be 
estimated as 



Fiis,t,u) 



s > 0,t > O. u > 0, 
{s + t + u = r) 



Fi(0,r,0)+ FiiO,t.u)+ Y Fiis,t,0)+ Y 

t>0,«>l s>l,t>0 min{s,u}>l 

O(Fi(0,r- 1,1)) = (^((^^"o-nogw), or 0(w""+"''-'-), 



(121) 



for rib < Ua or ni, > Ua, respectively, as ^ 0. Note that from Eq. pUlf) for = this estimation is valid for r = 
too. For r = nil + 2, it is estimated as 



C'(w"''-ilogw) (nb>2) 
0((log^)2) (nt^l) 



or 0(tj""^i), 



(122) 
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for ni, < Ua or Ub > ria, respectively, as — > 0. The summation of the second summand in Eq. (|117|l for < r < n^ + S 
is also estimated as 

r r 

X'lT F2{s,t,u)^ ^^2(s,i,0) = O(i^2(0,r,0)) = OK^+3-.)^ (^23) 

s > O. t > 0,u > 0, {s+t=r) 
{s + t + u = r) 

both for Tifc < Ha and for n?, > tIq, as w — > 0. The summation of the last summand in Eq. pi7f) for < r < + 1 is 
estimated as 

r r 

W Y F3{s,t,0) = Y F3{s,t,0)^O{F3{l,r~l,Qj) 

s >0,t>0,u>0, (s+t^r) 
{s -\- t -\- u — r) 

= 0(cj2"''-Mogw), or 0(0;""+"'-'-), (124) 
for Tib < Ha or Ua < Ub, respectively, as w — > 0. For r = + 2, the estimation is replaced by 

for Ub < Ua or na < rib, respectively, as — > 0. Then, by summarizing the above- noted estimations from Eqs. (|121|) 
to p25f) . and by taking into account the restriction (|96|l again, the proof of the lemma is completed. D 



B. The exceptional case of the first kind 



In this case, we first remember that from the discussion around Eq. (|5()|l it necessarily holds that Fi ^ 0, i.e., 
rifa = 1 in Eqs. (jHSJ- 

Lemma V.3 ; Assume that is an exceptional point of the first kind for H . Then the 0-th and the first derivative 
of R{z) are approximated by those of a finite series 



1 



A^zlogz 
that is, it is shown that 



A4z(logz)- 



+ A^QiAiQi + X''mQiriQi){QiriQi)-\ (126) 



R{z) 
1 



(QiAQi)-! 

A^z log z 

{QiFiQiy^Qi + A^QiAiOi + A27rzQiAOi)(girigi)-i 



A4z(logz)2 

= 0{z-^ {log z)-^) for r = 0, or 0{z-'^ {log z^^) for r = 1, 

as z — > 0. 

Proof : Let us first consider the quantity that 



dz"" 



R{z)- 
1 



1 



-(QiAQi)-! 



< 



A4z(logz)2 
d 



A^zlogz 

{QiriQi)-\Qi + X^QiAiQi + A27rzgiAQi)(QiAQi)"' 



dzJ 



R{z) - QiR{z)Qi 
1 



dz"- 



QiR{z)Qi 



}?z log z 



(QiAQi)-! 



A**z(logz)^ 



{QirxQiY^{Qi + A^QiAiQi + A27riQiriQi)(QiAgi)-' 



(127) 



(128) 
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For r = 0, the first term on the rhs of the above is estimated as follows: 

Riz) - QiR{z)Qij < \\QoR{z)Qn\\ + \\QoR{z)Qi\\ + jQiR{z)Qo\\ = 0(1), 
as z ^ 0, where Eqs. H68|) - (|70|l are used. For r = 1, one obtains 



dR{z) 


dQiR{z)Qi 


< 


dQoRiz)Qo 


+ 


dQoR{z)Qi 


+ 


dQiR{z)QQ 


dz 


dz 


dz 




dz 




dz 



(129) 



(130) 



In fact, by using expression H68|l the first term on the rhs of the above is estimated as follows: 
dQaR{z)Q„ 



= -QoRiz)Qo 

xQoRiz)Qo 
= O(logz). 



—z rfT^" - ^ 



(131) 

(132) 



Four derivatives in Eq. (|131|l have the same order, which can be shown from the use of Eqs. H()l|) . Ht)3|) . (|65|) . and 
(inH): We here note that 



O(logz), — — = 0(logz), — — = 0(logz). 



dz 



dz 



dz 



and from Eqs. Ht)l|) and H63|l 

^ = En' \^^^Qi{z + >^^[A{z)- log zr{z) + mr{z)]}Q,^E^,' 

^ 1^ f >?rc?^(^) , dr(z) driz) 

= 0(z-^(logz)-i). 

In the same way, the second term on the rhs of Eq. H13()() is also estimated as follows: 
Qo ( )Qi ^ _ I — JKL^^^ ^ £^00^ 1 °^ ) QiRQi + Ei^qEqiQiRQi 



(133) 



(134) 

(135) 
(136) 



\ c?z 



QiRQi 



= O(z-i), 

where we used Eq. (I69|l and the fact that 



-O(logz), — — = 0(logz). 



dz 



dz 



In a similar manner, we can also show 



dQiRiz)Qo 
dz 



O(z-i). 



(137) 
(138) 

(139) 
(140) 



Let us next consider the last term in Eq. H128|l . For r = 0, it reads 

1 



QiRiz)Qi 



X^z log z 



(QiAQi)-! 



1 



T^T^ ^{QiriQi)-\Qi + A^QiAiQi + AViOiAQi)(giAgi)"' 

A*z(logz)"^ 



< 



1 



-1 

" A^zlogz 

1 



(QiAOi)"' 



< 



X^z log z 



{QiriQi)-\Qi + X'^QiAiQi + A^^iQiAQOlQiAQi)"' 

-0-^00^-^01 ) ^ 

(QiAQi)"' 



A4z(logz)2 
-||(<3ii ^ £'i/£'io-Boo^-^oi)""^£']^/ - -E]^/ 

^ll(^) ^ 



1 



A^z log z 

i=2 ^ 

= 0(z-i(logz)-3), 
where in the second inequality we used Eq. H()2|) and that 



{QiriQi)-\Qi + X'^QiAiQi + A2^^QlAgl)(Ql^lgl)-^ 

-0(1) 



A4z(logz)2 

oo 



\E-^_-^^ EwEf^f^^ Eal\\\\E-^_-^^\ 



= o(i), 



1 ^ 11-^11 EiqEqq Eoi\\ 

as z — > 0. Substituting Eqs. H129() and (|144() into Eq. H128|) . we can obtain the estimation H127() for 
For r = 1, we can obtain 



dz 



QiR{z)Qi 



1 



A^zlogz 



(QiAQi)-' 



A4z(logz) 



2 (QiAgi)"'(Qi + A^gi^liQi + AViQiAQi)(QiAQi)"' 



< 



d_ 

dz 



1 



p-i 

" A2zlogz 



(QiAQi)"' 



1 



r(0iA0i)"H0i + A^giAiQi + A27rigiAQi)(QiAQi)"' 



< 



A4z(logz)2 

dz^^*^" ^ £^r/£^io-B(7o^£'oi)"^-Bfi^ - 
^" {^Qi + ^'[^(^) + ^«-r(z)]} Qij 



A4z2(logz) 



2 (QiAQi)"'(Qi + A^Qi^iQi + A2^iQiAQi)(0iAQi)~' 



-i?n |^A2(logz)Qir(z)Qi|i?rii - 



1 



-(QiAOi)"' 



1 



dz A4z2(logz)2 



dz A2z log z 

(giAQi)"'(gi + X^QiAiQi + AViQiAgi)(giAOi)" 



{(gii — -Eji^E'io-Eoo^-^oi) ^^'^ii}— 



{Qii - E^^ EiqEqq Eqi) E^^ 



d_ 

dz 

0(z-2(logz)-3), 
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where we used the expression for dE:[^/d z in Eq. (|134|l . Actually, the first term in Eq. (|148|l is estimated as 



A4z2(logz)2 



{QiriQi)-\Qi + A^QiAiQi + A27riQiAQi)(QiAgi)^' 



< 



11 



A^zlogz 



(QiAQi 



^Qi{^ + A2[A(z)+7rir(z)]}gi 



ll(giAQi)-i|| 



+ 



A2|zlogz| 
||(0iAQi)-i||2 



— Qi {z + \^[A{z) + 7rir(z)]} Qi 



1 



- T^n — (QiAOi)-' 
A^^zlogz 



A4z2(logz)2 

= 0(z-2(logz)-3), 



— Qi {z + X^[A{z) + ^^^(z)]} Qi - Qi(l + A^Ai + X^Tnri)Qi 



(150) 
(151) 
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as z 0. On the other hand, the second term in Eq. (|148|l is shghtly comphcatcd to evaluate: 

log 2)^ J 



< 



dz A4z2(log 

^(g.AQ.)-{^A^(iog.)Q.r(.)g,}3^(Q,AQ0 



-z 
d 1 
dz A^zlogz 
1 



(QiAQi)-! 



A^zlogz 
1 



(QiAQi)"'<{ ^A2(logz)Qir(z)QiUii ^ 



A"' 2; log z Idz J A'^zlogz 



A^z log z 
1 



< 



(A'^z log z)-^ \dz J 

(QiAQi)-' -A2(logz)Qir(z)gi — (giAQi)-' 

A^zlogz l_az J A^zlogz 

^iiT^l (QiAgi)-' -A2(logz)Qir(z)Qi V£;ii-^- (QiAQi)-' 

A^zlogz [dz J A^zlogz 

E^nT^l (QiAOi)-' -A2(logz)Qir(z)Qi 

A^zlogz 1^ dz J 

dzA2zlogz(«^^^^^)"| 

^ (giAQi)-'(:^A2(logz)Qir(z)QiUnT^4 (QiAOi)"' 

logz 1^ dz J A^zlogz 



J=2 
"A2^ 

d 1 



X^2 



+z 



+z 



1 



(A2zlogz)3 \d. 

1 



(A^zlogz)^ \dz 

+0(z-2(logz)-3) 

0(z-2(logz)-3). 



A^zlogz (QiriQi)-iz(Oi + A^QiAiQi + A^vriQi AQi)((3iAQi) 



(QiAQi)"' 



A'ziogz (girigi)-^(gi + A^gi^iOi + A^7rigiAgi)(giAgi) 



iiT^n — (QiAgi)-'(:^A2(iogz)gir(z)gij— i — 

A-^zlogz 1^ dz J A-^zlogz 



(152) 



(153) 



(154) 
(155) 
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In fact, the first term in Eq. H154|l reads 

(QiA0i)-'(;^A2 log zQ,r{z)Q A— ^ (giAQi)-' 

A^zlogz [dz J A-^zlogz 



< 



d 1 

o?0 X^z \ogz 



(QiAQi 



A2(logz)Qir(z)gi \ + — A2z(logz)Qi 



IKQiAQi 



z 



(A^z log z)2 

-(QiA0i)-'Qi^^0i + 0i 



(log z) 



x||(QiriQi)-i|| 



1 



A^zlogz 
= 0((zlogz)-i). 

The second term in Eq. (|154|l also reads 



1 



A^z log z 



(QiAQi)-!^ — A2(logz)Qir(z)Qi 



1 



^11 



A^zloez 



(QiAQi)-' 



1 



(A^zlogz)-^ \ dz 

< iKQiAQi)"' 



— A^zlogz {Qir,Qi)-\Qi + A^QiAiQi + X''7TtQiriQi)(QiriQi)-^ 



^A^(logz)Qir(z)Q] 



;(giAOi)~U-^Qi + A'Qi{A(z) - logz[r(z) - zA] + Tnr{z)}Qi} 

IKQiAgO^^II 



. dz 

< !l(QiAQi)-'|| 



A^zlogz (Qi + A'QiAiQi + A^TriQiAQi) 



(A^zlogz)^ 



— A^(iogz)Qir(z)gi f>(Qirigi) 



X \\-{zQi + A2Qi{A(z) - logz[r(z) - zA] + ^ir(z)}Qi} 
+z(Qi + A^giAiQi + A^TTzQiAQi)" ll^^i^i^i)"'" 



(QiAQi)-!] 



X ||z(Qi + X-'QiAiQi + X'niQiriQi) 
0((zlogz)-i), 



(A^zlogz)^ 
|^A2(logz)0ir(z)Qi|(QiAQi)-' 

(giAQi)-^!! 



^A^zlogz )Qi 



(A^zlogz)'^ 



(156) 

(157) 
(158) 



(159) 



(160) 
(161) 



as z ^ 0. The third term gives the same contribution to the order as the second one does. Furthermore, the last term 
in Eq. (|154|l comes from the estimations of the second, third, and last terms in Eq. p53|) . where each contributes the 
same order as 0(z^^(log z)^"^). Therefore, Eq. p55|) is proved. 
On the other hand, the third term in Eq. (|148|l reads 



{{Qii — E-^-^EiqEq^Eoi) ^ — Qii} 



dE 



11 



< 



\Eii Ei^Ef^^^ Eqi I 



1 ^ ll^ii^'^io-2'oo^-^oi I 



dz 



0{z-'), 



(162) 



as z ^ 0, where Eqs. H63() . H65() . H67|l . and 1136(1 are used. In the same way, the last term in Eq. (| 148(1 reads 
-E^n^^(Qii - -Bi/£'io£'oo^Ai)" 



< WE 



+ 



11 



e: 



{Qii - E^j^ EiqEq,^ Eqi) ^ 



^w^oo ^01 



.idE_ 



" dz 
O(z-Mogz), 



10 p-lp 
-C/QO ^01 



p-if? T^-i dEpi 
E,, E,oE,, 



dF~^ 
^11 -^10 — : -C/n 



dz 



(163) 
(164) 
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as z — > 0. By substituting Eqs. (|151|) . H155|) . (|162|l . and (|164|l into Eq. H148|l . one finally obtains Eq. H149|) . Wc can 

now show Eq. l(T77|) for r = 1 by setting Eqs. itHUI) and ffT^ into Eq. l(T^ . □ 

If we start with expression H72|) , we obtain the following lemma instead of Lemma IV.3I 



Lemma V.4 ; Assume that is an exceptional point of the first kind for H. Then the 0-th and the first derivative 
of R{z) are approximated by those of a finite series 



1 



A^z(log z — 27ri) 
1 



\'^z{\ogz- 2ni) 



(QiAQi)"' 

(QiriQi)-i(Qi + X^QiAiQi - \^mQiriQi){QiriQi)-\ 



that is, it is shown that 
d 



dz"" 



R{z) 



1 



A^z(log z — 2Tri) 
1 



M n ^^(QiAQi)"'(Qi + A^QiAiQi -A^7riQiriQi)(QiAQi)" 

A^z(logz — 27rz)^ 

0{z-^ {log z ~2Triy^) for r = Q, or 0(z"2 (log z - 27ri)"^) for r = 1, 



(165) 



(166) 



as z 0. 



VI. THE REDUCED TIME EVOLUTION OPERATOR 



In this section, we show that the reduced time evolution operator is expressed by the Fourier transform of the 
imaginary part of the reduced resolvent both in the regular case and the exceptional case of the first kind. We here 
define the reduced time evolution operator by the N x N matrix U{t) of the components Umn{t) ■= {m\Pe~'''*^ P\n), 
where P ~ £'((0,oo)) and {E{B)\B G B} is the spectral measure of H, which is a family of the projection operator. 
B is the Borel field of R. 

Lemma VI.l ; We assume that Eq. \2'^) holds so that there is no positive eigenvalue. Then, for the system with 
the rational form-factor 0), it holds that 

U{t) = - e-'*'^lmR+{uj)duj = lim - / e-'*'^Imi?+(cj)dw, (167) 



both in the regular case and the exceptional case of the first kind, where 



lmR+{uj) ^[R+{uj) - R-{uj)], (168) 
2i 

which is sometimes called the spectral density. 

Proof : Let us remember that the matrix U{t) is expressed by the spectral measure as 

Un.u{t) = I e-''^d{m\E{X)\n) = f e-"^d^„„,(A), (169) 

where E{B) is the matrix of the components {m\E{B)\n) . Therefore, what we first should do is to clarify the relation 
between Im_R+((jj) and E{X). Resorting to Stone's formula between E{B) and R{z), we clearly see 

hE{[a,b])+E{{a,b))]^ Urn [ [R{u; + le) - R{u; - ie)]du; , (170) 
2 t~,+Q 2m ,L 



for a, & S R with a < b. Under the assumption (|23|l . Lebesgue's dominated convergence theorem and the proof of 
Lemma IlIL II tell us that the exchange between the limit and the integration in Eq. (|170|l is allowed for [a, b] C (0, oo). 
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If [a,b] C (— oo, 0)\cr(J?), then = -R(w), and thus E{[a,b]) = E{{a,b)) = 0. In addition, by the continuity of 

ii^iuj), Eq. ((T7n|l tells us that E{{a}) = for all a > 0, which leads to 

1 /"^ 

E{{a, b)) = E{[a, b]) = - / Inii?+(a;)dw, (171) 

J a 

for all a, b with 6 > a > 0. 

Let us now consider the regular case and in particular the validity of the expression (|171|) for the interval including 
the origin. In this case, i?^(0) := linii^^+o -R^(<^) exists to be finite. Furthermore, limi^^ao R^{i^) = from Lemma 
IIII.2I Thus, is uniformly continuous on (0,(X)). Therefore, we can take the limit of Eq. (|171|l as a ^ +0 

to obtain liuia^+o E{[a,b]) = E{{0,b]). We next see that all components of lmR^{uj) are integrable, i.e., belong to 
i^((0, cx))). Suppose that |?/') G C^, then {ip\E{{0, X))\ip) is positive and a monotonically increasing function of A, 
and it is also differentiable in this case. Thus Eq. (|171|l tells us that {ip\lmR'^ {uj)\ip) > 0. In addition, 

I ~ 1 f°° ~ 

ll^ll^ > lim (V'|£^((0,A))|V') = lim - / {^j\lmR+ {uj)\i')duj = - {^j\lmR+ {uj)\i')duj. (172) 

Hence, from the monotonic convergence theorem, we see that {ip\lmR'^ {Lu)\ip) E i^((0, oo)). From this fact and the 
use of the polarization identity, we can prove that all components of the matrix Imi?+((x') are integrable. Thus, 
extending the rhs of Eq. H171(l to arbitrary B G {B G M\B C (0, oo)}, /g Imi?+(cj)du; defines a measure. We can now 
see from Eq. H171(l and from E. Hopf's extension theorem that 

E{B) = - / lmR+{uj)duj (173) 
TT Jb 

holds for all B E {B E M\B C (0,oo)}. Note that this expression means that the restriction of E„in{B) to {B E 
B|_B C (0, oo)} is absolutely continuous. Therefore, rewriting of U{t) in Eq. (|169|l into H167|) is straightforward. 

In the exceptional case of the first kind, from the assumption (|23|l . R^{uj) is continuous on (0, oo), while R^{uj) = 
0((u; logw)"^) as w ^ +0, so that it is not integrable around 0. See, Eq. H71(l . However, Imi?+(ijj) is of the order 
0(ij;^^(logijj)^^) from Lemmas IV. 31 and IV. 41 and thus it is integrable around 0. Hence, Eq. H173() holds again, and 
Eq. (|167|l is valid for this exceptional case. D 

We remark that in the case of no negative eigenvalues (point spectrum) of the restriction of l}{t) to the 

continuous energy spectrum is removed because in such a case P = I the identity. Furthermore, the connection 
between C/(i) and the observables is easily found, e.g., \{'ilj\U{t)\ij})\^ /\\P\tp)\\'^ for jV') e (or © {0}) is the 
survival probability of P\4') which is the probability of finding the system in the state P\i^) at the later time where 
P\i!) is just the decaying component of the initial state 



VII. THE ASYMPTOTIC EXPANSION OF THE REDUCED TIME EVOLUTION OPERATOR 

We can finally show the asymptotic formula for U{t) at long times for the rational form factors satisfying our 
assumptions. In the following, we assume that Eq. (|23|l holds, i.e., there is no positive eigenvalue. However, this is 
not explicitly mentioned in the statements of the theorems. Let us first consider the regular case. For this purpose, 
according to Lemma |V.2I we introduce the remainder F{uj) in the following way: 

-lTnR+{u) = X^iKm-^u^'^ir,,^ + ujr„^+i + oj^r„^+2){K{Q))-^ + F{uj), (174) 

TT 

for w > 0. 

Theorem VII. 1 ; Assume that is a regular point for H. For a system with the rational form factor ^ 
characterized by the positive integers Ha and that satisfy that ni, > 2 and Ua = 1, the reduced time evolution 
operator U{t) behaves asymptotically as 

u{t) - x^^^^iKm-TnAimr' + o{t--^-% (175) 

as t 00. When ni, — 1 and Ua > 1, the error term is replaced by 0(t~^ logt). 
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Proof : Wc first summarize the several properties of Imi?+(a;). By Lemma rV.2l we see that the remainder F{uj) in 
Eq. (|174|l is arbitrary-times differentiable. Particularly it holds that lim^^Q F{uj) /duj^ = for all r < rif,, and 



for nf, = 1 and > 1, or rif, > 2 and = 1, respectively, and 



O(logc^), or 0(1), 



^ ^' 1 0(loga;) (rib = 2) 



(176) 



(177) 



for rib = \ and > 1. or rif, > 2 and ria = 1, respectively, as w ^ 0. On the other hand, we see from Lemma llII.2l 
that {d/ dLoYlTa.R^ {lo) — 0{uj~'^~'^) as lo ^ oo. In particular, if m > 1, {d/ dioY'^lmR^ {uj) is integrable on [5, oo) for 
an arbitrary S > 0. 

Let us now split the integral in Eq. 11()7I) into two parts by writing 



Imi?+(w) 



{uj)lmR+ (lo) + (1 - 0(w))lm^+(w), 



(178) 



where (p G C^([0,oo)) and satisfies (j){u!) = 1 in a neighborhood of = 0. Such a function is realized by f{i^) ~ 
1 — Jq g{x)dx, where g{x) = h{x)/ J^h{x)dx and h{x) = exp(— l/[a-^ ~ (2; — d)'^]) {\x ~ d\ < a) or {\x — d\ > a) 
with d > a > 0. 

From Lemma 10.1 in Ref. ^4ind the above-mentioned discussion, we see that (1 — 0(w))lmi?+((jj) has a contribution 
of 0{t~'^) to U{t) for an arbitrary m > 1, i.e., this term decays faster than any negative power of t. 

On the other hand, the contribution of 0(tt')lmi?"'"(w) to U{t) gives the main part of the asymptotic expansion. 
Then, the coefficient of Fn^, J^nj^+i, and Fn^^2 is given by the form^ 



d^uj^ioj) 



dui^ 



(179) 



for all iV > 1 4- q, where q takes the value rib, + 1, or n?, -I- 2. We here used that 
d'^LoUitu 



duj^ 



min{fc,g} 

E 

J=0 





d^uj'i 


d''-^(j){uj) 






duj^ 







= r(i + (7)4g, 



(180) 



where r(l + n) = x"e ^dx is the gamma function. In addition, the remainder Rwit) is bounded above by 



\RNit)\ < 



°° d^L0l(t){uj) 



duj 



N 



e-"^'duj 



o{t--). 



(181) 



Note that since all derivatives of (j){uj) vanish in the neighborhood of w = 0, Eq. p79|) is valid for all > g + 1 and 
thus i?Ar(t) decays faster than any negative power of t. Furthermore, we understand, by applying Eq. (|B1(I in Lemma 
B.ll directlv to 4>{uj)F{uj) with the discussion in the first part of this section, that the contribution of the Fourier 
transform of the remainder (j){uj)F{uj) to U (t) is 



0{t-^\ogt) or Oit-"^"-^), 



(182) 



for rift = 1 and tIq > 1, or > 2 and ria = 1, respectively, as ti; — > 0, where we used the formula of the indefinite 
integral that J {[ogioYdio = a;[(logix>)^ — 2\oguj + 2]. Summarizing the above-noted results, we finally obtain that 



< 



r(l + nh) 



F„. 



(l-(/)(cj))Imi?+(w)e-**"^dw 
O(t-Mogt) or 0(i-"''-2), 



r(2 


+ rib) 




"6+2 




-/ 




Jo 



r(3 + nfc) 

^"6 + 1 + /..N„(,+3 ■'"(.+2 



{Lu)F{u;)e-'"^duj 



(X(0))-i 

+ oit-^) 



(183) 



for nf, = 1 and ria > 1, or rib > 2 and = 1, respectively, as i — s- 00, where 0{t ^) is due to the contribution from 
i?Ar(i). This is just the asymptotic expansion of U{t) in the statement. D 
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It is worth noting that if we resort to Lemma FV. II instead of Eqs. (I176|l and (|177|l . we have 

for nil ~ I and tIq > 1, or Tif, > 2 and Ha ~ I, respectively. However, in the latter case, we see that the Fourier 
transform of (I){uj)F{uj) gives the contribution of the order 0{t~'^''~^), which is just the same order as that coming 
from the dominant one. Hence, we can only obtain an useless estimation. 

We next show the asymptotic formula for U (t) at long times for a system with an exceptional point of the first 
kind. To this end, we write 1t[iR'^{lu) with the remainder F{lu) again as follows: 

ilmi?+(^) = -— ^ -iQiriQ.r' + Fiu). (185) 

TT A^CL;(ioga;)^ 

Theorem VII. 2 ; Assume that is an exceptional point of the first kind for H, which necessarily imposes that 
nij — 1. Then, the reduced time evolution operator U(t) for the rational form factor 0) behaves asymptotically as 

Uit) = ^^^(giAQi)-' + 0((logt)-2), (186) 

as t — !■ cx). 

Proof : Let us first look over the some properties of Imi?+(a;) again. By Lemmas IV.3I and IV. 41 we see that the 
remainder F{uj) in Eq. H185|l is arbitrary-times differentiable, satisfies that F{uj) = 0(a;^^(logcj)^^), and dF{uj)/du! = 
0(aj~^(logw)~'^), as a; ^ +0. On the other hand, we see from Lemma flll.2l that {d/dLjjyimR~^{uj) = 0(uj~'^~^) as 
Lu — > oo. In particular, if m > 1, {d/ dLo)"^lmR'^ (uj) is integrable on [S, oo) for an arbitrary S > 0. 

We now split the integral in Eq. H167|l into two parts as in Eq. H178|) again using the C^^-function 4>{oj). From 
Lemma 10.1 in Ref. ^^nd the discussion mentioned above, (1 — (f>{uj))lmR'^ (uj) has a contribution of 0{t~™) to U{t) 
for an arbitrary m > 1. On the other hand, the contribution of (/)(ti;)Imi?^((jj) to U(t) gives the dominant part of the 
asymptotic expansion. Then, the dominant time dependence of the asymptote of U{t) follows from Lemma [6.21 that 



(^(a;)(w(loga;)2)-ie-"" duj = (logt)-^ + 0{{logty^). (187) 







Furthermore, the contribution of the Fourier transform of the remainder (f)(u})F{ijj) to U{t) can be estimated by the 
similar manner to Lemma [6.21 rather than Lemma FB. II By setting cr{uj) ~ i^(a;)e~**'^ instead in the proof of Lemma 
.21 we can apply it to this case, and we have 

F{uj)(j){uj)e-''^du = - Yim {ia){u) + {-lf [ {1^ (j){u)^-p^dLo. (188) 



uj^W Jo duj^ 



Then, corresponding to Eq. ljB10|l . we have 



F{lo - ir])e-*''dr] 







+ E{t)<C Uuj-ir])-^[log{uj~iT])]-^\e-'^'idT] + E{t), (189) 







with an appropriate constant C. Note that in this procedure, i?+(w) is analytically continued to the lower plane of 
the second Riemann sheet, while R~{uj) still remains in the lower plane of the first Riemann sheet. Then, both are 
ensured to contribute the remainders of the same order to F{uj — irj) in the above integral from Lemmas IV. 31 and 
IV. 41 The remainder term E{t) in Eq. (|189() that gives the order of 0{e~^*) for some 7 > is responsible for the 
possible poles of R'^{u}) continued to the second Riemann sheet, the number of which are guaranteed to be finite from 
the analytic Fredholm theorem^k and Lemma fill. 21 for the continued i?+((jj). Thus, it follows from Lemma FB . 21 that 
limij_»+o ll(-^f)('^)ll = 0((logt)~^). By the same argument as in Lemma fB. 21 one also sees that the remainder term 
in Eq. H188|l is of the order of 0(t~^~^^). Summarizing these arguments, we can finish the proof of the theorem. D 
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VIII. CONCLUDING REMARKS 



We have rigorously derived the asymptotic formula of the reduced time evolution operator for the A'^-level Friedrichs 
model in the context of the zero energy resonancsi^ both for the regular case and the exceptional case of the first kind. 
Then, in the latter case, the logarithmically slow decay proportional to {\ogt)~^ has been found, and the expansion 
coefficient has been explicitly presented by the projection operator associated with the zero energy eigenstates of the 
total Hamiltonian, which is an extended state not belonging to the Hilbert space. We note that the decay involving the 
logarithmic function expressed by t~^{logt)'' [j = 1,2,... and fc = 0, ±1, . . .) can occur in the short range potential 
systems in the even dimensional spacei^ It should be noted that a realization of the exceptional cases require the 
parameters, e.g., the coupling constant A, to take such special values that the matrix K{0) in Eq. H44|l has a zero 
eigenvalue. In addition, some of the form factors Wn(w) have to behave as |'(;„(tj)p ~ c„w around lo = 0. In other 
words, if all of them behave as |w„([j)p ~ c„a;*" with q„ > 2, the exceptional case of the first kind never occurs though 
that of the second kind could happen. These circumstances explain how the exceptional cases arc surely exceptional. 
The presented results also enable us to calculate the asymptotic formula for the survival probability of an arbitrary 
initial state 1-0) localized over the N discrete levels. If we choose the special initial state to satisfy rn^{K{{)))~^\il}) = 
in Eq. H175|) or = in Eq. H186|) . our estimations are useless and other decay laws could appear iS£ The 

long time behavior of the reduced time evolution operator for the exceptional case of the second and the third kind 
are not examined. As is expected, in the former case, the non decaying component associated with the localized zero 
energy eigenstate wih appear due to the divergent behavior of Q2R{z)Q2 = 0{z~^) in Eq. The latter case can 

occur in the A^-level cases of the model only for iV > 3, which yields a more complicated situation. In the whole of the 
paper, we assumed that there is no bound eigenstate with a positive eigenencrgy. This situation is actually realized in 
the weak coupling cases^^*^ However, its compatibility with the existence of the extended zero energy eigenstate is 
still not clear in the multilevel cases (except the single level case). The emergence of the logarithmic decay (logt)^^ 
is just due to the logarithmic energy dependence of the self energy S{uj) and it comes from the assumption ^ where 
|t;„(tj)P ^ Cnio'^"' with a positive integer qn is required. Therefore, if we choose another type of form factor, it is not 
necessary for such a slow decay to occur even in the exceptional case^i*iSii^ However, we stress that our assumption 
is often satisfied by actual systemsii^ii^ The experimental realization of the exceptional case requires the setup of 
parameters like ~ A-^A, where A is a typical cutoff constant. This seems in a strong coupling region to be naturally 
satisfiedii^, and hence it could be suggested to invoke the artificial quantum structures for a realization. 
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APPENDIX A: CHARACTERISTICS OF SELF ENERGY FOR THE RATIONAL FORM FACTOR 

Lemma A.l Suppose that "qiw) is a rational function, i.e., it is expressed by ri{uj) = t:{uj)/p{lu), where n(uj) and 
p{lo) are the polynomials of the degree m and n, respectively. Furthermore, we assume that n > m + 1 and p{z) has 
no zeros in [0,oo). Then 



for all C G C\([0,oo) U {ak}^^i), where Uk is a pole of ruk-th order of rj{z), N is the number of such poles, and Pn{C) 
is a polynomial of C of the degree not greater than n . For ( = |C|e*^ with < 6 < 27r, we define —C, = ICk**^ with 

— TT < < TT. 

Proof : From the fundamental theorems for the complex functions, it holds that 

!M_du.^~''^Res(^Aogi-z),z^aX (A2) 
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where aN+i = C- Then the residue at z = for k < N is deduced to exphcitly 



{mk - 1) 



2 mk~'- 

~Ty. ^ 



^TOfc - 1 



ri{z) 



For z = Cj which is a simple pole, the residue becomes 

viz) 



Res 



l0g(-z),Z = C =7?(C)l0g(-C). 



(A3) 



(A4) 



Therefore, by setting Eqs. (|A3(I and (|A4() into Eq. (|A2|I . one obtains Eq. ljAl|l . and the proof is completed. D 

Lemma A. 2 Suppose that the function rj^uj) belonging to L^{[0,oo)) is of the form 

7]{uj) := ujPr{uj), (A5) 

where p > 1 and r(uj) is a C^-function defined in [0, oo). Then it holds that both ri{uj)/uj and r]{uj)/uj'^ £ L^{[0,oo)), 
and 



lim — 



duj 





-i 


Jo ^ 





ri{uj) 



duj. 



(A6) 



Proof : From the proof of Proposition 3.2.2 in Ref. 10 , the principal value of the integral on the rhs is written by 
the absolutely integrable function as follows 



P 



LU- E^ 



duj 



ri{uj) - i']{E)(ps{uj - E) 
oj-E 



dio. 



(A7) 



for all E > 0, where ips{uj) is a C^-function with support [—(5,(5] {0 < 6 < E), even with respect to the origin, and 
such that ^ps{0) = 1. In the following, we choose such that (ps{u}) = exp[l — 1/(1 — (uj/S)'^)] for lu £ (—(5,(5) or 
otherwise, and S = E/2. Then, the proof of Eq. I|A6|I is equivalent to that of 



lim 



7]{uj) r]{uj) - r]iE)ips{uJ - E) 



ri{uj) 



duj = 0, 



LU — E 

which will be shown in the following. We note that the above-mentioned integrand can be rewritten as 

ri{uj) i]{E)ips{uj - E) 



-E 



lu^{lu-E) E{lu-E) 
r,{E)ipsiuj - E) Eriiiu) - ujri{E)>fs{^ - E) 



Euo 



uj^{uj-E) 



(A8) 

(A9) 
(AlO) 



We also put h = {0,E/2], h = {E/2,2,E/2), and h = [iE /2,oo). 

Let us first consider the case that e /i U /s. Then, since ^psi^^ — E) = 0, we can use Eq. I|A9|I to estimate the 
integrand in Eq. IjASp . which reads 



E 



r]{uj) 



lo^{lo~E) 



< 2 


77(w) 




W2 



(All) 



where the rhs is absolutely integrable and independent of E. Furthermore, it follows that 

lim£;^+o £'x/iu/3(w)77(u;)/[w(a; - E)] = for every uj e (0, oo), where Xhuhi^) = 1 (t^ e /i U I3) or 
(lu g I2), being the characteristic function. Thus, by the dominated convergence theorem, we can see that 



lim 



E- 



ri{uj) 



..... , - 0. T^\(^^ = 



EXhuhi^) 



r]{uj) 



lu^{lu-E) 



du) = 0. 



(A12) 
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Next, for w € /2, we can use Eq. ljA10|) . The integration of the first term of Eq. (|A10|I is estimated by 



I2 



7jiE)ips{LU - E) 

Elu 



duj 



< 



Ey2 



(ps{u! — E)duj 



I2 



E 



ipi {x)dx —* 0, 



(A13) 



as i? ^ +0, because ~ 0((jjP) where p > 1. The second term of Eq. HA10() is also estimated with the 

decomposition 

\ET^{iu) - u;t]{E)M^ -E)\<\{E- u;)r^{iu)\ + \lu{7^{u;) - TjiE))\ + - - E)\. (A14) 

The integral corresponding the first term on the rhs of the above is evaluated as 

\{E - u;)7^{u;)\ 



LU'^IUJ^EI 



-duo 0, 



(A15) 



as E' ^ +0, because of the fact 'q{uj)/uj'^ e i^([0,oo)). The integral corresponding the second term is also evaluated 
as 



uj^\lo-E\ 



div < (In 3) sup \t]'{uj)\ 0, 



(A16) 



as -E +0; because of the assumption on ri{uj), where is the derivative of r](uj). The integral corresponding the 

last term on the rhs of Eq. (|A14() is also estimated as 



J2 Lu-'luj-El 6 |^|<i 



a,s E ^ +0. Thus, we see from Eqs. (|JT5)l . l|JT6jl . and lPl7)l . 

lim / Eriiiu) - ujiiE)^s{uJ - E) 



duj = 0. 



uJ^{uj-E) 

Equations HA12|I . (|A13() . and (| A 18(1 mean the completion of the proof of Eq. (|A8() . □ 



(A17) 



(A18) 



APPENDIX B: ASYMPTOTIC EXPANSION OF THE FOURIER INTEGRALS 

We have to estimate the integrals of the form U{t) = e^^*^F{X)d\ in which F{\) = identically either for small 
A > or for large A where F is supposed to take values in an arbitrary Banach space B. The following lemma is 
essentially the same as Lemma 10.2 in Ref. 1^ 

Lemma B.l .• Suppose that F{X) ^ for X > a > 0, (F <E C''+^{S, oo;B)), fC'+i) e L^{d,oo;B) for any 6 > 
and for an integer k > 0, and that F^^^O) = for j < k — 1. Then 

\\Um<^{ WF^'HmdX+^J'^ sup \\F^'^+^Hf^)\\dx] , (Bl) 



for all t > 7r/a. Here F^'^^X) denotes the k-th derivative of F(X) and so forth. 

Proof : By extending F by F{X) = to A < 0, we obtain a function F on (—00,00) with F^'^^ e i^(— 00, 00; B). 
Then we have that 

||E(^)(A + /i)-FW(A)||(iA = ir + / ) ||E(^)(A + /i) -FW(A)||(iA (B2) 



r2h pa rX-\-h 

<2 ||FW(A)|ldA+ / dA / ||F(*'-+i)(^)|ld/i (B3) 



2h 

< 2 / \\F'^^\X)\\dX + h I sup \\F^''+^\^i)\\dX. (B4) 

Jo Jh fj.e[\.X+h] 

By noting that e^**^ = _g-»t(A-7r/t)^ g^gg ^-^^^ ^j^g jj-^g Qf |jB2() is just an upper bound of the Fourier transform 
of 2F^''\ Remember that the Fourier transform of F^*^' is equal to {it)^U{t) under the assumption of the lemma, i.e., 
F'^-')(0) = for J < — 1, then the desired result follows. D 
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Lemma B.2 Suppose that 4> € C|5"([0,oo)) and satisfies 4>{uj) ^ 1 in a neighborhood of to = 0. It then holds that 
for any positive integer q > 2 and N > 1, 



+0((logt)i-''-^) (B5) 

as t — 5- cxD . 

Proof : We first put cr(ijj) = ijj~^(loga;)~*e~**'^ and introduce the indefinite integral operator / as2S 

s-i(logs)-«e-'*"(is, (B6) 
wliere c is an arbitrary complex number. We can also recursively show 

(/V)(w) = — ^— / (w- s)'=-is-i(logs)-«e-''*"ds. (B7) 



Then, repeating the partial integration, we obtain 

POO 

/ uj-\\oguj)-muj)e-'"^duj = ~{ia){ij)cb{uj)\^^o + (B8) 
Jo 

for all iV > 1 with 

RN{t) = (-1)^ / {i^a){uj)d^4>{oj)/doj^duj, (B9) 



where we used the fact that d^ (t){Ld) / du)'' — at oj = for any fc > 0, and d^ (piuj) / duj^ = at = cx) for any fc > 0. 
We now choose c = w — ioo and change the variable as s := w — iij, which leads to 



{ia){uj) = ie-'"^ / {lo - M7)-^[log(w - iT])]-''e''"dr]. (BIO) 
Jo 

Then, we can use the dominated convergence theorem to obtain 

\im {i<j){oj) = i (-i?7)-i[log(-i?7)]-?e-*''d77. (Bll) 

Here we used the fact that there is a positive number < l/(e''V2) such that 

\i.-^,r[^o,i.-^,)^^e-^^\ < { ' i^^;!/^""' ^\^'<:;i'^ 

[x,<^Mv-'\iogir' + C]e~'\ (B12) 



< 



for all < uj < loq and all < ry < oo, where x-n<'^oi''l) ^ 1 V < or otherwise, and C is an appropriate 
constant. The existence of such a C is ensured by the fact that log(a; — irj) has no zeros in the rectangular region 
{uj — iij\0 < UJ < ujoi^o < ^ < oo}, and its modulus diverges as r/ oo. The function on the rhs of Eq. I|B12|I is 
integrable, so that the use of the dominated convergence theorem is valid. To evaluate the asymptotic behavior of 
Eq. (|B11|) . putting e = (logt)~^ and trj = ^, one obtains 

hm = ^{-ey / -[1 - e\og{-^0]'-'' e^^d^ (B13) 

uj^+O l-q \d^ ) 



/ [l-6logH0]^-^e-«d^ 



(B14) 



1 — o I ^ \ j I \dv 2 
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where we used the formulas that /(e) = Y.^S^ f^^\0) / j\ + [t^ /{N - - u)'^-^ f^'^\eu)du with /(e) = 

[1 - elog(-iO]^"'', and x''-'^e-^'^{\ogxy dx = d'[ii-''T{v)]l dv^ . The last integral in Eq. (|bT5)| turns out to be 
finite because wc have 



{l-u) 



N-l 



Jo [Jo [l-e^xlog(-zO]''+^ 
where we used that 



du 



[log(-zO]^e-«dC 



< 



[(7r/2)^ + (logO 
(7r/2)9+^ 



,f:\2]q/2+N 



e-^d^, 



|1 - eulog(-ze)P = (1 - eulogO' + (W2)^ > (^,^2)^ + (logCr ' 
Let us now evaluate the upperbound of <t){lu)\. Using the estimation (|B12|I . we have 



|(/^a)(a;)| 



1 



< 



(iV-1)! 
1 

W^Jo 



\if-^\\ogujo\-'' + Cr;^-i]e'*"d7? = 0{t-^+'). 



(B16) 
(B17) 

(B18) 
(B19) 



Substituting this result into the error term Riy{t) in Eq. (|B9|I . we see that RN{t) = 0{t ^+i) for any integer N > 2. 
This proves the statement of the lemma. D 
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